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ABSTRACT 

Context. To determine the physical parameters of a transiting planet and its host star from photometric and spectroscopic analysis, it 
is essential to independently measure the stellar mass. This is often achieved by the use of evolutionary tracks and isochrones, but the 
mass result is only as reliable as the models used. 

Aims. The recent paper by Torres et al (2009) showed that accurate values for stellar masses and radii could be obtained from a 
calibration using T^ff, log g and [Fe/H]. We investigate whether a similarly good calibration can be obtained by substituting log p - 
the fundamental parameter measured for the host star of a transiting planet - for log g, and apply this to star-exoplanet systems. 
Methods. We perform a polynomial fit to stellar binary data provided in Torres et al (2009) to obtain the stellar mass and radius 
as functions of Tgff, log p and [Fe/H], with uncertainties on the fit produced from a Monte Carlo analysis. We apply the resulting 
equations to measurements for seventeen SuperWASP host stars, and also demonstrate the application of the calibration in a Markov 
Chain Monte Carlo analysis to obtain accurate system parameters where spectroscopic estimates of effective stellar temperature and 
metallicity are available. 

Results. We show that the calibration using log p produces accurate values for the stellar masses and radii; we obtain masses and radii 
of the SuperWASP stars in good agreement with isochrone analysis results. We ascertain that the mass calibration is robust against 
uncertainties resulting from poor photometry, although a good estimate of stellar radius requires good-quality transit light curve to 
determine the duration of ingress and egress. 
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1. Introduction 

There are currently over 400 known exoplanets, of which more 
than 60 transit their host star^ This important transiting sub- 
set are the only planets for which the orbital inclination, and 
hence precise stellar and planetary parameters, may be deter- 
mined. The fundamental parameters found for the host star and 
transiting planet are stellar density (see below) and planetary 
surface gravity ( Southworth et al. 2004 1. To convert these into 
values for the radii of both, it is necessary to find the stellar 
mass. This is often aiTived at iteratively via deriving a stellar 
density from the lightcurve analysis and a stellar effective tem- 
perature from spectroscopy and using these with model evolu- 
tionary tracks and isochrones of appropriate metallicity to find 
a stellar mass and age (Soz zetti et al.|2007) l. Further photomet- 
ric and spectroscopic analysis may then be performed to aiTive 
at final values for the masses and radii of the star and planet, 
see e.g.lHebb et al. (2009 1. The resulting values for masses and 
radii are therefore only as reliable as the evolutionary models 
used. A recent study by Southworth (2009) highlighted the fact 
that discrepancies between different sets of evolutionary models 
represent the dominant source of systematic uncertainty in plan- 
etary parameters. For example, they find that the spread of mass 
values obtained for HD 209458 using different models is around 
4%. 

Here we develop a new one-step approach to determining 
the masses of exoplanet host stars from their effective tempera- 



tures, metallicities and photometric bulk densities. We base our 
method on the recent study by |Torres et al.| ( [2009[ ) of the masses 
and radii of a large sample of well-characterised main-sequence 
stars belonging to non-interacting, eclipsing spectroscopic bina- 
ries. ToiTes et al. (2009) showed that accurate stellar masses and 
radii could be obtained using a calibration of stellar surface grav- 
ity, effective temperature and metallicity. They used a set of well- 
determined measurements of log g, Tgff, [Fe/H], M and R from 
binary stars to obtain coefficients that allow mass and radius to 
be calculated directly for any normal star, without isochrone fit- 
ting. 

Recently, the use of log p in place of log g in the determina- 
tion of star-planet system parameters has become widespread, 

( |2008] l. 



see for example Sozzetti et al. (20071, Winn et al. 



ISozzetti et al. ( 2009 ) and [ Fernandez et al.| ( |2009| l. Where high 
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quality photometric data can be obtained of the transit event, the 
stellar parameters can be obtained more precisely using the stel- 
lar density value derived from the lightcurve than using the stel- 
lar surface gravity value from spectral analysis (Sozzetti et al. 
2007|). 

In Section 2 we review the methodology for determining 
exoplanet host-star densities from the transit geometry. We re- 
determine the mass and radius calibrations of |Torres et al.| ( |2009l l 
using their data, and obtain comparably tight mass and radius 
calibrations using log p in place of log g. In Section 3 we ap- 
ply the method to the host stars of several transiting planets for 
which isochrone mass determinations have been published re- 
cently. In Section 4 we show how the method can be incorpo- 
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rated directly in a Markov-chain Monte Carlo (MCMC) analysis, 
to give the stellar mass as a derived parameter. 



2. Analysis 



Table 1. Coefficients for Mass and Radius fits. 



We used the tabulated data of Torres et al. (2009i to perform 
our calibration. Those data consist of 19 binary systems, i.e. 38 
stars, for which the metallicity is known, after excluding systems 
that contain pre-Main Sequence stars. We fit a similar polyno- 
mial calibration, replacing log g with log p, and apply it to star- 
exoplanet systems. The stellar density can be obtained directly 
from only photometric measurements via fitting of a transit event 
( Seager & Mallen-Ornelas 2003 1: the ratio of semi-major axis to 
stellar radius depends on the ratio of transit duration to orbital 
period via 



a 



(1 + ^FKFf-bHi - , 



^2nT 



1/2 



(1) 



where AF = {F - F,ra„sii)/F, b is the impact parameter, T is the 
transit duration and P is the orbital period. 
Combining this with Kepler's Third Law, 



2„i 



G{M, + Mp) 



(2) 



where G is the gravitational constant, leads to an expression for 
the stellar density ( [Seager & Mallen-Ornelas|2003) 



Pg 



4k^ 



P^G 



(1 + yfKFf -b^{\- sin^'f) 



■ill 



(3) 



since M, >> Mp (this is strictly true for planets on circular or- 
bits only). Thus the stellar density may be obtained directly from 
parameters measurable from a high-quality lightcurve: the dura- 
tion and depth of transit, the impact parameter and the orbital 
period. 

We used a Singular Value Decomposition (SVD) fit, 
weighted by error on the mass or radius measurements, to ob- 
tain coefficients on significant variables. The final fit gives coef- 
ficients for X {- log{Teff) - 4.1), X^, log p, log p^, log p^ and 
[Fe/H] for mass, and X, log p and [Fe/H] for radius (the sec- 
ondary terms are insignificant for the radius fit). Thus the mass 
or radius may be computed by 

logM = al+a2X+a3X^+a4logp+ailogp'^+a(,\ogp^+a^[Fe/H](4■) 



logR = bi+ b2X + bj,logp + b4[Fe/H] 



(5) 



To obtain Icr errors on the coefficients, we carried out a 
Monte Carlo analysis of 50,000 runs in which, for each run, each 
value for T^ff, log p, [Fe/H] and M or R were perturbed ran- 
domly on a gaussian with standard deviation the Icr error for 
that measurement. The set of 50,000 resulting coefficients were 
recorded, and the standard deviation of each coefficient within 
that set gave the error value for that coefficient. The fitted and 
error values for each of the mass and radius fits are given in 
Table [T] The resulting scatter in fitted less measured values is 
log M = 0.023 and log R = 0.009. Figure [T] shows calibrated 
versus measured mass and radius. 



Mass, a, 



Radius, bj 



const 


0.458 ± 


0.017 


0.150 


± 0.002 


X 


1.430 ± 


0.019 


0.434 


± 0.005 




0.329 ± 


0.128 






logp 


-0.042 ± 


0.021 


-0.381 


± 0.002 


log p~ 


0.067 ± 


0.019 






log 


0.010 ± 


0.004 






[Fe/H] 


0.044 ± 


0.019 


0.012 


+ 0.004 



3. Application to WASP host stars 

Seventeen WASP host stars have now been published. We use 
values for Tgjf and metallicity obtained from spectral analy- 
ses with the coefficients given in Table [T] in Equations |4] and |5] 
to calculate the mass and radius of each, and compare these to 
the values obtained from isochrone interpolations. WASP-1 data 
was obtained from from [Cameron et al. (2007), [Stempels et aT 



( [2007| , [Charbonneau et al.[ pOOTTand .Shporer et al.| ([20071l, 



WASP-2 data from [Cameron et al. ( 2()07|l an d" Charbonneau [ 
[etaL[ ( |2007] l WASP-3 d ata from| PollaccoetaL[ p008l), WASP- 
4 data from [Wilson et "al . (2008), WASP-5 data from Anderson] 
^taTjlWES, WASP-6 d ata from Gill on et al.[ ( |2009l ), WASP-7 
data from Hel her et al.l (f20 09b), WASP- 10 data from |Johnson[ 
etH] ([200911 a nd Christian eTaTI (|2009]l WASP- 11 data from 



West et al 



2009b^ and |Bakos et a 



p009]l, WASP- 12 data from 



Hebb et al.| ^20091 , WASP- 13 da t a from[SEIIen et al.[ ( |2009] l, 
W ASP- 14 data from [Joshi et al.[ ( |2009] l, WASP- 15 da t a from 
iWest et al.[ ( [2009a| l , WASP- 16 data from [Lister et al.[ ( [2009] l, 
WASP- 17 data from[Anderson et al.[ ( |2009l l, WASP- 18 data from 
IHe llieret al.[ ( [2009a 1) and WASP-19 data from'Hebb e t aL] ( |20T0l l. 
Table [2[ shows that the agreement in results from the two meth- 
ods is very good, demonstrating that computing stellar masses 
and radii from these coefficients is a valid, and simple, alterna- 
tive to isochrone interpolation. 

The only values in Table [2] that do not quite agree, within 
errors, are the results for the radius of WASP- 10. [Johnson et al.] 



(20091 find a radius of 0.70^^ ;^ 7?©, whereas the calibration pro- 



duces R = 0.60'^qI\Ro. WASP-10 is an unusu al host star, with 
a high density of 3.10po Johnson et al. [[200911 and a high level 
of activity ( [Smith e t al. 2009 ). Such calibration discrepancies in 



low-mass, high-activity stars are discussed in Section 6 o fjTorres 
|etaL] ( |2059l l 



4. Application to Markov Chain Monte Carlo 
Analysis 

Having established that our modification of the Torres calibra- 
tion yields stellar masses that agree well with those obtained 
from isochrone fitting, we now describe its implementation in 
our Markov-chain Monte Carlo (MCMC) parameter fitting code. 
This is followed by a discussion of a Main Sequence mass-radius 
constraint generally imposed on the parameter fitting, and ob- 
jects for which this constraint is removed. We also discuss the 
effect on the final mass and radius values of uncertainty in ec- 
centricity due to imperfect radial velocity and photometry. 

The MCMC analysis, described in [Collier Cameron et al.[ 
(2007 1 and Pollacco et al. ( 20081 1, derives star-planet system pa- 
rameters from simultaneous modelling of stellar lightcurve and 
radial velocity measurements. The host star's radial velocity mo- 
tion is parametrised by the radial velocity amplitude, the centre- 
of-mass velocity, the orbital eccentricity and the longitude of 
periastron, while the transit is modelled using the formulation 
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(a) (b) 

Fig. 2. Shows the scatter in isochrone (a) mass and (b) radius values versus values from MCMC analysis for 17 SuperWASP host 
stars. 



of Mandel & Agol' ('2002') incorporating the 4-coefficient limb- 
darkening model of Claret (2000 ). The form of the transit profile 
is determined by five parameters: the epoch of transit, orbital 
period, duration and depth of transit and the impact parameter 
of orbit. The best model is found using a constrained optimiza- 
tion of for both photometric and radial velocity data, com- 
bined with Bayesian priors on epoch, period, transit duration 
and depth, impact parameter, stellar mass, radial velocity semi- 
amplitude, eccentricity and the longitude of periastron. 

Here, we adapted this analysis to take the spectral T^ff and 
metallicity as input values when modelling all available pho- 
tometry and radial velocity observations of a host star. The 
MCMC code determines the stellar density value for the calcu- 
lations from the observations at each step in the chain, now us- 



ing Bayesian priors on the temperature and metallicity also, and 
uses the equations and coefficients above to calculate the stellar 
mass. These modifications make the MCMC analysis more ro- 
bust since the stellar mass is now a derived quantity in the chain 
instead of a jump parameter constrained by a prior 

Where there is high-quality follow-up photometry of a transit 
event, the duration of ingress and egress, and hence the impact 
parameter, are tightly constrained. However, where such pho- 
tometry is not available and thus the duration of ingress and 
egress cannot be accurately measured, an additional constraint 
is needed in the MCMC analysis. In such a case, the transit 
ingress and egress durations are overestimated in modelling the 
photometry since the sharp transitions to ingress and egress be- 
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Table 2. Comparison of mass and radius values obtained from isochrone fitting with those from the calibrated equations and from 
an MCMC analysis using those equations. 















Isochrone 




Fitted 




MCMC 




ID 


p. 


Teff 




Metallicity 


M, 


R. 


M, 


R. 


M. 


R. 




Po 


K 








Me 


Re 


Me 


Re 


Me 


Re 


WASP-1 


0.38+!;!!! 

-u.Ul 


6110 


±45 


0.26 


±0.03 


J 2j+om 

-0.07 


1.49+!!!1^ 

-0.04 


2 21+0-05 

-0.07 


1.50+!!-"? 


1.20+"!13 

-O.Oj 


1.41+!!-!' 

-U.14 


WASP-2 


1.51+" l« 

— U.13 


5150 


±80 


-0.08 


±0.08 


0.82+I!'i 

— U.U/ 


0.82+"!!« 


0.88+!!-°^ 

— U.Uo 


0.82+°-"3 

— U.Ui 


0.88+""2 

— U.Uz 


0.89+"-°^ 

— U.U4 


WASP-3 


-0.05 


6400 


± 100 


0.00 


±0.20 


J 24+0 08 

-0.08 


1 31+"'2 

-0.05 


2 21+0.05 

-0.07 


1 32+**-*^^ 

-0.09 


1.22+""^ 

-0.04 


1^30+"-°i 

-U.U4 


WASP-4 


1.09+!! 

-0.09 


5410 


±240 


0.00 


±0.20 


0.90+!:"' 

-0.0/ 


1 15+"-'* 

-0.28 


0.95+!!-"* 

-O.Uo 


0.95+"-"^ 

-U.03 


0.91+"!!' 

-0.05 


0.88+!!-"2 

-U.UZ 


WASP-5 


1 .08^n H 

-u.II 


5700 


± 150 


0.00 


±0.20 


o.99+!:"^ 

-0.08 


0.97+""* 

-O.Oo 


i.or!!-"? 

-0.0/ 


0.97+"-"' 

-U.09 


1 ■01^04 

-0.04 


i.oo+"-"5 

-U.UZ 


WASP-6 


J 34+0.11 

—0. 10 


5450 


± 100 


-0.20 


±0.09 


0.83+!!"' 

— u.uy 


o.85+!!!!3 

— U.U4 


0.93+!lI!^ 

—U.Uo 


0.87+!!-"2 

— U.UJ 


0.93+""^ 

—U.UZ 


0.93+"-"^ 

— U.U4 


WASP-7 


o.6r?,f, 


6400 


± 100 


0.00 


±0.10 


1 9^+0.04 

— U.UJS 


1.23+!!"' 

-0. lo 


1.19+"!!^ 

—U.U/ 


1.23+I!-"^ 

— U. 12 


i.2o+"!!3 

—U.UJ 


1 21+0-04 
—0.04 


WASP- 10 


3.10+!!!!' 

-0.09 


4675 


± 100 


0.03 


±0.20 


0.75+!!"t 

-0.0 J 


o.7o+!!!ll 

-0.01 


0.80+"!!^ 

-U.U / 


0.60+!!-"! 

-U.Ol 


0.79+""? 

-U.UJ 


0.69+"!!' 

-U.UJ 


WASP- 11 


0.69+!:?? 

— u. 1 1 


4800 


± 100 


0.00 


±0.20 


0.82+!!"3 

—U.U J 


0.75+""2 

— U.UZ 


0.84+!!-°^ 

—U.Uo 


0.76+"-"2 

—U.UJ 


o.8o+"!!3 

—U.UZ 


0.79+!!-°3 

—U.UZ 


WASP- 12 


o.35+!:°3 

— U.UJ 


6300+^"" 

"■^""-100 


30+" "^ 

"■-^"-0.15 


1 33^*^-^^ 

-0.05 


1 57+0.07 

—0.07 


2 27+0-06 

—0.08 


1 C7+0.05 
—U.U/ 


1.28+"!!' 

-U.U5 


1.60+!!;!! 

— U. lU 




Q 43+0.12 

"■^-'-0.10 


5826 


± 100 


0.00 


±0.20 


1 03+" " 

•"-^-0.09 


1 34+"" 


1 10+"-°* 

^■^"-0.08 


1 40+"- 

^■^"-0.17 


1 11 +""3 

^■^^-0.03 


1 94+0.02 

^■^^-0.02 


WASP- 14 


"■-'^-0.06 


6475 


± 100 


0.00 


±0.20 


1 21+" " 


1 Q 1 +0.07 
■-^ -0.07 


1 9^+0.05 
■ ■-^-0.07 


1 34+0.05 

■- ^ -0.07 


1 T 1 +0.05 
^■■-^^-0.10 


1 4Q+0.07 
^■^'-0.29 


WASP- 15 


n Q7+0.04 
"■-"-0.04 


6300 


± 100 


-0.17 


±0.11 


1 IQ+O.IO 

^•^'^-o.io 


1 48+""' 

^■^°-0.07 


1 20+"-"5 

'^■'^"-0.07 


1 c 9 +0.05 
— ^ -0.07 


1 1 Q+0.03 


1 Q7+0.10 
i--"_0.14 


WASP- 16 


1 21+"" 


5700 


± 150 


0.01 


±0.10 


1 00+" "5 

^ -""-0.07 


94+0.04 

"■'^-0.05 


1 or"-"' 


93+" "' 

"■'-^-0.05 


1 00+" "3 

^■""-0.03 


97+0 03 

"■"-0.03 


WASP- 17 


67+" '* 

"■"'-0.13 


6550 


± 100 


-0.25 


±0.09 


1 16+"'^ 

^•^"-0.12 


1 20+""'* 

^■^"-0.08 


1 1Q+0.05 
^■^'-0.07 


1 9 Q +0.09 
^■^■-'-0.11 


I 1 q+0.03 

^■^=^-0.03 


1 1 7+0.07 
^ ■ ^ ' -0.06 


WASP- 18 


71+""* 

"■ ' ^-0.10 


6400 


± 100 


0.00 


±0.09 


1 94+0.04 

^•^^-0.04 


1 20+" "* 

^■^"-0.10 


1 lQ+0.05 
^■^'-0.06 


1 1 9+0.06 
^■^^-0.06 


1 99+0.03 
^■^ -0.03 


1 97+0.04 
^■^'-0.04 


WASP- 19 


1 13+"'2 

^■^-'-0.12 


5500 


± 100 


0.02 


±0.09 


95+" "9 

"•'-'-0.10 


94+" 0' 

"■'^-0.07 


n Q7+0.04 
"■="-0 06 


94+0 03 
"■'^-0.04 


93+0 02 

"■'-^-0.02 


n Q 1 +0.03 
U.81_(,_o2 



come blurred, allowing a shallower slope to be fitted. This leads 
to an underestimation of the impact parameter, b, from 



1 - ^^6- 



(6) 



from Winn| ( |2009 1 where 6 is the transit depth, T is the total 
transit duration and t is the partial, flat transit duration. This 
leads to an overestimation of the stellar radius, Rt from 



acosi 



(7) 



for a circular orbit, modified from |Win"n| ( |2009[ l, where a is semi- 
major axis and ; is orbital inclination. To avoid this overestima- 
tion of stellar radius, the MCMC-fitted stellar radii are generally 
constrained to reasonably closely follow a Main Sequence re- 
lationship to the steUar mass, i.e. R^ - (Cox 2000). To 
achieve this, a Bayesian prior is imposed on the stellar radius 
within the MCMC analysis such that a penalty is added to a 
chain step with a stellar radius different to that expected for the 
stellar mass value of that step, i.e. 



O.S\2 



2 

Xadd 



(Rsten - K l 



(8) 



where Rstep is the stellar radius value for the current step in the 
Monte Carlo chain, M» is the stellar mass of this chain step and 
cr^ is the estimated uncertainty on the power-law estimate of the 
stellar radius. This penalty is added to the value for the 
chain step, giving the model a lower chance of being accepted. 

However, WASP-1, 12 and 15 have this Main Sequence con- 
straint relaxed due to being more evolved stars, each with an es- 
timated age greater than their expected Main Sequence lifetime. 



WASP-1 and 12 are late F-type stars with an estimated main se- 
quence lifetime of < IGyr but estimated ages of around 2 Gyr. 
The high-quality photometry of Charbonneau et al. I^^QffTj and 
IShporeret al.lpOTj i (for WASP-1) and |Hebb et al. 



(2009 I (for 



WASP- 12) constrains the impact parameter and stel ar density 
satisfactorily, allowing a good estimate of stellar age, and at the 
same time removing the need for the additional MCMC con- 
straint. WASP- 15 is an F5 star of around 3.9 Gyr compared to 
a normal F5 age of 3.6 Gyr, also with high-quality photometry 
( West et al.|2009a i which again constrains the impact parameter 
and stellar density. Restraining the radii of these stars to lower 
values with the Main Sequence constraint on leads to an over- 
estimation of their density. In addition, we used the very high- 
quality photometry for WASP- 10 given in Johnson et al. ( 2009| l, 
and therefore had no need for the Main Sequence constraint in 
that case. For all other objects, without extremely high quality 
follow-up data and expected to be on the Main Sequence, the 
constraint was kept. 

With high precision photometry the precise shape of the tran- 
sit would help constrain the eccentricity of the planetary orbit 
through the measured duration of transit. A planet moving on an 
eccentric orbit has its velocity modified by ( l+esm(L))/( Vl - e^) 
compared to an identical planet on a circular orbit. Therefore for 
a planet on an eccentric orbit, the duration of transit, T, given 
by the length of the chord crossed in front of the star divided 
by the planet's velocity, is modified by ( Vl - e^)/(l + esinw) 
compared to the duration in a circular orbit. Without very high 
quality photometry, a source of uncertainty in the stellar density 
can arise from the fitting of the eccentricity of the planet's orbit 
from the radial velocity measurements. With stellar density as an 
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Table 3. Result of fixing e = for WASP- 13 



e 


b 


p. 


M. 




0-18oT 
0.0 


15° °2 
09°'° 


90° 
0.58°;°2 


1 05°°^^ 

1 11°°^^ 


1 05°°' 
1 940.02 

^■^^0.02 



input, this leads to uncertainty in the results for stellar mass and 
radius using the calibrated equations and coefficients. This may 
especially be an issue for low mass planets which have a lower 
ratio of radial velocity amplitude to scatter than more massive 
planets, making an accurate determination of orbital eccentricity 
more difficult. Using the simplified version of the equation for 
density. Equation 19 from Seager & Mallen-Omelas ( |2003[ l, us- 
ing T-njP « 1 so that sinx * x, and then including this factor 
gives 



Pg 



32^ 
Gn 



1 + esinw\ 



3/2 



(9) 



where t is the partial transit duration consisting of the 'flat' part 
only. 

As an example, WASP- 13b is a planet of mass 0.46My or- 
biting a Gl type star with a r adial velocity semi-amplitude 
of 56OTi ' (Skillen et al. 2009 |l. The only transit photometry 



available is from Super WASP and the 0.87m James Gregory 
Telescope in St Andrews, Scotland. Running the MCMC anal- 
ysis as described above results in an eccentricity value of e = 
0.18 + 0.05, with;^'^ = 15136. However, holding the eccentric- 
ity fixed at and repeating the analysis leads to an output with 
a very slightly higher value of 15140, so clearly the larger 
eccentricity value is not a secure result. The difference in result- 
ing parameters are given in Table [3j showing that there is only a 
small effect of about six percent on the final stellar mass value 
despite a change of over 50% in the stellar density. Therefore it 
may be concluded that the effect of the uncertainty in eccentric- 
ity on the final fitted mass value is small. The stellar radius is 
altered by around 18%, due to the change in density value. 

The results of the MCMC analysis on each of the 17 
SuperWASP host stars are also presented in Table |2] and are 
shown in Figure |2] Eccentricities are held fixed at in all cases 
where this was done in the original papers WASP-3, 4, 5, 7, 1 1, 
13, 15, 16 and 19). Almost all mass and radius values from the 
MCMC analysis agree with the isochrone values, within errors; 
the MCMC mass and radius ranges of WASP-6 and the radius 
range of WASP- 19 do not quite overlap with the isochrone val- 
ues, but all these values agree at the 2cr level. 



the stellar density estimate which in turn requires good knowlege 
of the impact parameter. Thus in establishing planet radii there 
is no substitute for good quality photometry, though the Main 
Sequence prior can provide a useful additional constraint if the 
star can be shown via independent means to be unevolved. 
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5. Summary 

We have presented a new calibration for stellar masses and radii 
based on stellar effective temperature, metallicity and stellar 
density. We have shown that the resulting equations provide a 
good fit to data for 38 stars from [Torres et al.| ( |2009| l, and also 
to values for masses and radii of exoplanet host stars obtained 
from isochrone analyses. We have demonstrated that accurate 
stellar masses may be obtained for such exoplanet host stars via 
a Markov-chain Monte Carlo analysis of photometric and spec- 
troscopic data, using spectroscopically determined temperatures 
and metallicities as input. 

Even where poor photometry yields an uncertain estimate of 
stellar density, the mass estimate from the calibration is encour- 
agingly robust. However, the stellar radius depends strongly on 



